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A simultaneous multi-point, multi-component velocimeter using interferometric 
detection of the Doppler shift of Rayleigh, Mie, and Rayleigh-Brillouin scattered light in 
supersonic flow is described. The system uses up to three sets of collection optics and one 
beam combiner for the reference laser light to form a single collimated beam. The planar 
Fabry-Perot interferometer used in the imaging mode for frequency detection preserves the 
spatial distribution of the signal reasonably well. Single-pulse multi-points measurements of 
up to two orthogonal and one non-orthogonal components of velocity in a Mach 2 free jet 
were performed to demonstrate the technique. The average velocity measurements show a 
close agreement with the CFD calculations using the VULCAN code. 


I. Introduction 

T he gas flow velocity is an important parameter to measure in scientific, engineering and industrial applications. 

In aerodynamics, the velocity field dominates effects such as heat transfer and mass transport. Novel and 
improved methods of measuring velocity are sought for gases because existing methods do not work in all 
environments and/or require undesirable modifications to the flow, such as the need to seed particles or gasses into 
the flow. It is desirable to have a measurement system that can be used to determine flow properties, particularly 
velocity, but also temperature and density, in unseeded gas flows. This method should ideally work at both high and 
low temperatures, at multiple points, should measure multiple velocity components, and should be instantaneous. In 
this paper we present a variation of the Rayleigh/Mie scattering method that meets these requirements. 

One approach to measure Rayleigh and Mie scattering spectra is to use atomic and molecular filters. In Doppler 
Global Velocimetry 1,2 (DGV) two cameras view the flow, one of which looks through an Iodine filter (I 2 ). If the 
laser frequency is kept constant, variations around this frequency due to Doppler shift correspond to variations in 
transmission. Using a calibration curve between transmissivity and velocity, up to three components of velocity can 
then be calculated. Although with the evolution of the seeding and computing techniques this method reached a high 
degree of accuracy, the principal limitation is that particles need to be seeded into flow. Only at low gas temperature 
can the DGV methodology be used in the Rayleigh scattering regime (scattering from molecules without seeding). 
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Filtered Rayleigh Scattering 2 3 (FRS) has been 
developed for unseeded measurement of velocity, 
temperature, and density in airflows and can be 
used for 2D mapping. The filter can also be used 
to suppress the background scattering. However, 
the I 2 cell has limitations in that the absorption 
line from the spectrum that is used for filtering 
has a small bandwidth. Consequently, this 
approach has a low velocity dynamic range and 
does not work at high temperatures (>1000 K) 
when the Rayleigh spectrum is much broader than 
that of the filter. 

A second approach permitting a larger 
dynamic range of measurable velocities and 
working at higher temperatures is the 
Interferometric Rayleigh Scattering method. This 
technique has previously been used for single- 
point velocity, density, and temperature Figure 1. Theoretical spectrum of Rayleigh scattered light 
measurements. Fabry-Perot and confocal etalons for different values of the y = p/rjKa parameter using the 
were used for Rayleigh signal analysis in Tend S6 model. Reprint from Reference 8 by permission 
conjunction with cooled CCD, intensified CCD, of the author and of the American Institute of Aeronautics 
or photo-multiplier array that image the and Astronautics, Inc. 
interferogram. 4,5,6 Online referencing using a 

planar Fabry-Perot increased the precision of this method. 7 Furthermore, it was shown that multiple points could be 
measured simultaneously with a Rayleigh scattering interferometric approach, but only one component of velocity 
was measured. 8 This approach was used in some measurements without in-line referencing of the laser light, 
decreasing the instrument precision. On the other hand, in order to increase the number of photons collected this past 
work has generally been limited to excitation with continuous-wave lasers and integration for long times, resulting 
in single-point time-averaged measurements of velocity. Three velocity components have been measured 
simultaneously by using combined instruments of different types (Rayleigh and LDV) and having different methods 
of obtaining the velocity. 9,10 This approach is very expensive, complicated and hard to use. 

It has also been demonstrated that two orthogonal velocity components can be measured simultaneously at one 
point using one interferometer by reflecting back the probing laser beam 1 1 or at multiple points by using a second 
optical cavity. 12 However, this approach results in directional ambiguity of the flow velocity vector. Prior 
knowledge of the approximate magnitude and sign of the velocity components can remove this measurement 
ambiguity. The method described in Reference 12 removes the ambiguity by slightly detuning the secondary cavity. 

This paper describes an optical system for gas velocity measurement using laser scattering from molecules or 
naturally-occurring particles analyzed for Doppler shift by a planar Fabry-Perot interferometer. Velocity 
measurements are performed simultaneously at multiple spatial locations. Multiple velocity components are 
measured. This system will be combined in the future with Coherent Anti-Stokes Raman Spectroscopy (CARS) 
thermometry 13 to allow simultaneous temperature, species concentration, and velocity measurements in supersonic 
combustion flows. 

The velocimetry system collects the Rayleigh scattered light from different angles. This light is mixed together 
with the unshifted laser light before it is passed through a planar Fabry-Perot interferometer for analysis. At the 
output of the interferometer, a high-sensitivity CCD camera images the interference fringe pattern. This pattern 
contains both the Rayleigh scattered spectrum from multiple collecting angles and the reference laser spectrum 
information. The shape of the scattering spectrum is function of the parameter y = a/Ka, where a = pip is the 
collision frequency (p is the gas pressure and p is the shear viscosity) and K a is the acoustic spatial frequency (K is 
the wave vector and a is the speed of sound). For low density gases, when the parameter y « 1, the Rayleigh 
spectrum can be approximated with a Gaussian shaped spectrum. The peak of the spectrum is Doppler shifted by a 
frequency proportional to the bulk velocity, the spectral width is proportional to the square root of the gas 
temperature, and its amplitude is proportional to the gas density. For higher density gases, a more detailed 
theoretical model is required to describe the Rayleigh scattering spectrum. As shown in Fig. 1, the Rayleigh 
spectrum for higher density gases, calculated using the Tenti S6 kinetic theory model, 14 includes three Lorentzian- 
shaped peaks in the spectral shape. To measure velocity, the Rayleigh spectrum is analyzed for Doppler shift 
applying a fitting technique and using the laser spectrum as reference. 
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Figure 2. Experimental setup for simultaneous measurement of three components of velocity. 


II. Instrument Configuration 
A. Collecting Optics and Scattering Geometry 

The optical setup of the instrument is illustrated in Fig. 2. The vertically polarized laser beam with the wave 
vector k 0 is focused in a small region at the measuring location. The Doppler shift frequency of the scattered light 
collected in the direction of k si (perpendicular to k 0 ) is given by Avi = (k s i - k 0 ) ■ V, where the k si is the wave 
vector of the scattered light. The direction of (k s i - k 0 ) vector gives the direction of the velocity component Vi being 
measured. If a second collecting optics is placed along the same optical axis, but on the opposite side of the probed 
region, the scattering wave vector is k s2 = - k sl , where |k sl | = |k s2 | = |k 0 |. A mirror M r is used to reflect this collected 
light along the optical axis, sending it back to the interferometer. This light shows a Doppler shift Av 2 = (k s2 - k 0 ) ■ 
V and the velocity component V 2 is measured along the vector (k sl + k 0 ), i.e. the two measured components are 
orthogonal. By slightly tuning the mirror M r in the vertical direction the image will contain two parallel lines 
situated symmetrically around the center of the interferogram. 12 It will be shown that, each horizontal line represents 
one measurement direction. 

The Rayleigh scattered light in the k sl direction is collected and collimated by the L 2 lens. The beam diameter is 
reduced to match the interferometer’s input-aperture using lenses L 3 and L 4 . In the k s2 direction, the light is collected 
and collimated by the lens Li and reflected back toward the interferometer by the mirror M r . The third component of 
velocity V 3 with the wave vector k s3 is measured in a vertical direction in a plane containing k sl and k s2 but 
perpendicular to the plane defined by the vectors k sl and k 0 . The angle between k si and k s3 is about 33°. The light is 
collected and collimated by the lens L 6 , and the lenses L 7 and L 8 reduce the beam diameter. After the half wave plate 
HP 2 rotates the beam polarization to horizontal, the mirror M 3 reflects the collimated beam toward the polarization 
dependent beam combiner BCi. A Dove prism (DP) is used to rotate the laser beam image with an angle that is 
convenient for image processing (usually 90°). It can be positioned either on the Vi (and V 2 ) or V 3 optical path. High 
quality achromatic lenses were used for collecting and collimating optics and ^/20 flatness mirrors for M r and M 3 . 
Thereafter, the light is mixed together with laser light (used as reference frequency) before it is passed through the 
Fabry-Perot interferometer (F-P) using the beam combiner BC 2 . 

At the output of the interferometer, a high-sensitivity CCD camera images the interference fringe pattern thus 
the probed region through the etalon. 8 Therefore, the interferogram contains spatial as well as spectral information. 
Each set of fringes along the laser path in the image corresponds to a different measurement location in the flow. 
Because it is required that the light be collimated and perpendicular to the etalon input mirror M b the field of view 
is restricted to a small region at the focusing point of lens Li. Dedicated imaging software processes each 
interferogram, extracting the Doppler shift of the Rayleigh/Mie signal relative to the reference laser frequency and 
determining the velocity values. An interferogram containing both the Rayleigh scattered signal (a straight vertical 
line pattern) and the reference laser frequency information (concentric rings) is shown in Fig. 4b. 
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The air-gap adjustable planar interferometer 
was set for a free spectral range (FSR) either 3 or 
5 GHz, depending on the velocity range to be 
measured. If the magnitude of the velocity of 
interest would have only positive or only negative 
values, as for example in the case of axial 
velocity, the instrument could use the full FSR for 
measurement without ambiguity. In turbulent 
regions, where the radial velocity can change 
sign, only half of the FSR can be used for 
measurement without ambiguity. 

Due to the fact that our image recording 
technique employed on-line referencing, it wasn’t 
necessary to use the interferometer temperature 
controller, neither was it required to take special 
care against vibrations. 7,12 


B. Pulse Stretcher 

An increase in the number of photons 
collected by the receiving optics is desirable since 
Rayleigh scattering signals are low. This can be 
achieved in two ways: by decreasing the laser 
wavelength or by increasing the laser energy. In 
our experiments we have chosen to increase the 
number of photons up to 4 times by increasing the 
light source energy to 300 mJ/pulse using a 
Nd:YAG laser at 20 Hz. The laser beam diameter 
was about 10 mm and we used an f/#4 lens with 
the focal length of 400 mm to focus the beam at 
the measurement location. The breakdown of the 
gas limited the energy of the laser pulse at about 
65 mJ/pulse for this configuration. To avoid this 
problem, the 10 nsec pulse of the investigation laser beam was stretched to more than 40 nsec. 

The optical arrangement of the pulse stretcher designed to obtain three consecutive laser pulses of almost 
identical amplitude is shown in Fig. 3. The system contains two beam splitters BSi (68%T) and BS 2 (38%T) and six 
identical flat mirrors M x _ 6 to create three delay lines of 0.3, 1.0, and 1.0 m respectively. The glass plates GPi and 
GP 2 are used to compensate for the change in the optical path by refraction of the laser beam that occurs at the beam 
splitters. A small portion of the vertically polarized light reflected from one of the compensating plates is collected 
by an optical fiber to be used as reference (ROF in Fig. 2). Two cross-polarizers (not shown) are used to adjust the 
reference intensity collected from a stray reflection from one of the compensating plates. At full input laser power, 
the pulse stretcher can deliver 210 mJ/pulse with an efficiency of about 70%. Just over half of this energy is 
delivered to the probing volume because of losses in the transmission optics. The output beam pulse temporal shape 
differs slightly from the design shape shown in Fig. 3b. The tolerance of the commercially available beam splitters, 
and the setup constraint on the path lengths of the delay lines contributed to this difference. 

C. Beam Combining Configurations 

Figure 4 shows four possible combinations of mixing laser light and Rayleigh scattered light collected from 
different angles. The circular fringe pattern shown in Fig. 4a is generated by the reference laser frequency with no 
signal present. The vertical line in Fig. 4b, the image of the laser beam is used when multipoint measurement of one 
component of velocity (Vi) is desired. The positioning of the Rayleigh scattered light is symmetric on the 
interferogram and allows fast processing by averaging a small number of rows (reference) and columns (Vi), which 
pass through the center of the interferogram. In Fig. 4c the signal from two orthogonal velocity components, Vi 
(left) and V 2 (right) are injected off axis in the vertical direction. A configuration as shown in Fig. 4d with Vi 
(vertical on axis) and V 3 (135°) can be used to measure two non-orthogonal velocity components. The advantage of 
Fig. 4d over 4c is potentially faster processing because the interferogram can be read directly along radii at angles of 
0°, 45°, 90° and 135°. Figure 4e shows a combination of two orthogonal velocity components in the horizontal 
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(a) 




BSi 




Figure 3. (a), optical setup of the pulse stretcher using 
three delay lines, and (b) the design shape of the output 
beam. 



(a) (b) (c) (d) (e) 

Figure 4. Combining on the same interferogram Rayleigh scattered light from different angles with 
reference laser light: (a), reference laser light; (b), multipoint one component of velocity; (c), and (d) 
multipoint two components of velocity; (e), multipoint three components of velocity. 


direction (Vi and V 2 ) and one non-orthogonal in the vertical direction (V 3 ). Since in the case of an interferogram of 
type (c), and (e) the interferogram must be processed entirely, the processing time is considerably higher than in (b) 
and (d). 

Although the intensity of the reference laser light must be monitored carefully to be of the order of the received 
scattered light, the system does not require special devices for filtering this radiation. As will be shown later, 
referencing the scattered light with unshifted laser light in a single interferogram reduces the velocity measurement 
uncertainty. 


D. Image Processing 

Image processing software is used to process 
and analyze the entire interferogram for the 
Doppler shift of the Rayleigh signal and the 
measured free spectral range (distance between 
consecutive ring patterns), and to solve for 
velocity. The image processing methodology 
consists of linearization of the concentric 
interference patterns from a planar Fabry-Perot 
interferogram and subsequent analysis. The 
important variables are the center of the 
interferogram, the free spectral range of the 
interferometer, and the maximum fringe shift of 
the Doppler shifted signal. An interferogram- 
processing example of a pattern containing the 
reference laser frequency information (concentric 
rings, as shown in Fig. 4a) is shown in Fig. 5. 
There are four steps: 1. The image center is found. 
2. The image is mapped by reading it in r-6 plane, 
where r is the radius from the center and 6 is the 
angle, using the original (r, 0) coordinate as the 
new (x, y) coordinate (the resulting image is 
shown in Fig. 5a). 3. The coordinates of the 
obtained data are rescaled in accordance to the 
linear pixel configuration of the CCD camera. For 
small reading angles 6 \ the conversion described 
in steps 2 and 3 can be simplified by averaging a 
low number of lines or columns symmetrically 
along the interferogram diameter. 4. The fringe 
spacing is transformed by the radius relation R p+1 
= (R2 + 2fi/k)l/2 between two consecutive rings 



r(a.u.) 



(b) r < au > 

Figure 5. Transformation of concentric interference 
patterns from a Fabry-Perot interferometer (Fig. 4a) to a) 
parallel lines, and b) parallel lines of equal spacing 
spectra. 
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1.2 




Figure 6. a) Fabry-Perot spectrum of reference laser light, and b) two linearized Rayleigh-Brillouin spectra 
of low Doppler shift. The spectra of Mie or Rayleigh scattering for y « 1 are similar with the Gaussian 
shape spectrum shown in (a). The RMS in Fig. b) shows the goodness of fit when forcing three Lorentzians 
and one Gaussian as the fit function. 


of the Fabry-Perot interferogram, where R is the radius of the ring, p is the fringe number, k is the order of 
interference, and / is the focal length of the fringe forming lens. This step completes the conversion of the image, 
namely stretching data over concentric circles into data distributed over parallel lines of equal length. The final 
image, Fig. 5b can be processed directly in 2D by autocorrelation, 6 7 or by binning the lines in a single average 
spectrum around a chosen angle 0 and applying a fitting technique. In the present work, a fitting technique using as 
a function of the flow conditions, multiple Gaussians for Rayleigh and Mie scattering 12 (Fig. 6a), and a combination 
of Lorentzians and Gaussians for Rayleigh-Brillouin scattering (Fig. 6b), was applied to the linearized spectrum to 
determine the free spectral range (FSR) and Doppler shift frequency Af for each fringe. The velocity component was 
then calculated knowing the experimental geometry and using the equation Af= (k s -k 0 )F7/l where, k is the incident 
laser wavelength, k 0 and k s are the incident and the scattered wave vectors, and V is the gas velocity. 


III. Velocity Measurement 


A. Mach 2 Cold Jet 

Experiments to demonstrate the technique in a noisy environment were carried out in a small-scale Mach 2 jet 
facility at NASA Langley Research Center. The jet was formed with a convergent-divergent nozzle designed by the 
method of characteristics to produce a supersonic flow of M = 2 for a y= 1 .3 (M is the Mach number of the flow and 
y is the perfect gas constant). As designed, it had a throat diameter of 7.51 mm and an exit diameter of 10 mm; it 
was subsequently found that, due to manufacturing error, the throat diameter was 7.66 mm (a fact significant in that 
the CFD work described below used the “as-designed” contour). The jet discharged into the atmosphere. No special 
care was taken for drying or filtering neither the air supply nor the ambient air. Naturally occurring Mie scattering 
from water clusters condensation was present only in the mixing layer region. Although the system can measure 
both Mie and Rayleigh scattering by adjusting the CCD camera gain, this need for adjustment made it difficult to 
obtain a continuous scan across the jet diameter. A camera with automatic gain control or with a larger dynamic 
range would address this problem. 

The jet was mounted on a 2D stage, allowing axial and radial scans of the jet at the maximum speed of 0.55 
mm/sec in increments of 0.1 mm. The instantaneous shadowgraph images of the flow (approximately 1 psec 
exposure time) taken at different stagnation chamber pressure and temperature values are presented in Fig. 7. Three 
experiments were performed with nozzle plenum pressures of 5.5, 8.6 and 13.6 atm and plenum temperature of 
about 285 K. This corresponds to an over-expanded case, a perfectly expanded case and an under-expanded case. 
We probed a volume of about 1.5 mm length and 0.2 mm diameter perpendicular to the flow direction. The spatial 
resolution of each measurement point was about 0.2 x 0.2 x 0.2 mm 3 . The minimum distance from the nozzle exit 
that could be measured was at about 0.5 mm from the exit plane. Since the interferogram contains spatial as well as 
spectral information, the separation distance between the first two consecutive fringe orders on one side of the 
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interferogram corresponds to about 0.2 mm in the 
flow and the diameter of the first order fringe to 
about 0.5 mm in the flow. 

For comparison, numerical simulations of the 
experimental flow field were performed using the 
VULCAN CFD code. This program solves the 
discretized Navier-Stokes equations using a 
structured-grid finite-volume formulation of the 
governing equations. 15 Inviscid fluxes were 
computed using Roe's upwind flux difference split 
scheme and the turbulence was modeled with the 
K-omega model of Wilcox. The gas was assumed 
to be calorically perfect with a ratio of specific 
heats of 1.4. Three calculations were performed 
using the pressures and temperatures in the 
plenum chamber identical to those from the 
experiment. The grid used in the calculations 
started in the plenum upstream of the nozzle 
throat and extended 61 nozzle-exit diameters 
downstream of the nozzle exit. In the radial 
direction the outer boundary was more than 40 
nozzle-exit diameters from the centerline. There 
were 609 grid points between the nozzle exit and 
the downstream outflow boundary, with 401 of 
those points in the first 5 nozzle-exit diameters 
downstream of the nozzle exit. This produced a 
streamwise grid spacing of 0.05 mm at the nozzle 
exit and a grid spacing of 0.2 mm at a 0.5 exit- 
diameters downstream of the nozzle exit. There 
were 801 points in the radial direction between 
the centerline and the outer boundary. At the 
nozzle exit, there were 161 evenly-spaced points 
between the centerline and one nozzle radius, 
giving a radial grid spacing of 0.03125 mm. In 
order to reduce numerical difficulties associated 
with solving a jet flowing into quiescent air, a 
slight pressure gradient was imposed across the 
external domain. This, along with the entrainment 
from the jet produced an external flow of 
approximately 7 m/sec. The radial (v) and axial 
calculations is presented in Fig. 8. The coordinate x 
jet with the origin situated axially in the nozzle exit 



76.2 


50.8 


25.4 




Figure 7. Instantaneous shadowgraph images of the jet at 
the stagnation pressure and temperature of (a) 5.5 atm 
286.5 K, (b) 8.6 atm and 286.9 K, and (c) 13.6 atm and 283 
K. The exposure time is about 1 ps. 
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Figure 8. Numerical simulations of the flow field using the 
VULCAN CFD code. The conditions are as for the over- 
expanded case shown in Figure 7a. 

(u) velocity distribution of the over-expanded jet using these 
is in the streamwise direction and y in the radial direction of the 
plane. 


B. Radial velocity along jet axis 

The interferogram configuration shown in Fig. 4b was chosen for this measurement. In one test the jet was 
scanned axially in the streamwise direction at constant speed while recording interferograms at the laser pulse 
frequency of 20 Hz for 90 seconds. To compute the radial flow velocity (Vi), five rows were averaged through the 
center of the interferograms. This data was converted from the spatial domain (pixel) to the frequency domain as 
described earlier and in more detail in reference [9]. One spectrum processed in this way is shown in Fig. 6b. A 
combined function of one Gaussian and three Lorentzians were fitted to this data to determine with sub-pixel 
resolution the location of the peaks and implicitly the frequency. The low amplitude peaks in each spectrum are the 
sidebands of the Rayleigh-Brillouin spectra (see also the theoretical spectrum of Fig. 1). The two peaks shown in 
Fig. 6b (corresponding to the left side of the interferogram) show the simultaneous measurement at two locations in 
the flow situated at about 0.2 mm apart. Similarly, there are two locations on the right side of the interferogram, up 
to a total of four points in each image. The spectrum of the reference frequency (laser frequency) used to identify the 
zero-velocity fringe location (extracted from a perpendicular direction in respect to the previous spectrum) was fitted 
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with Gaussian functions (Fig. 6a). A series of 
1000 interferograms (8000 spectra) were 
processed in this way. Fig. 9 shows the average of 
the radial velocity measured simultaneously at 
four points on the jet axis. The mean of the 
measured FSR was calculated using all data and 
its variation around the mean was converted to 
velocity as if it would be a Doppler shifted signal. 

The result was plotted with an offset of -300 
m/sec in the same figure for comparison. This 
fluctuation in the measured FSR was found to be 
less than 1% of the full range of the instrument. 

Translated in velocity, this corresponds to a 
standard deviation of about 25 m/sec. 

Because the measurement volume was off 
axis by about 0.5 mm, the radial velocity in the streamwise direction is not zero. The data shows the variation in 
velocity when the measurement volume passes through the shock structure. Considering that the effect of turbulence 
was minimal in the region investigated, the velocity spikes seen in the plot are mainly due to the errors introduced 
by the on-line fitting method. 

For comparison, an average of the radial velocity at four points of the measurement volume scan was extracted 
from the CFD data. The result was plotted against the measured data and is shown in the Fig. 9 as a red line. 
Although the measured data shows a change from maximum to zero velocity faster than the CFD data in the 
expansion region of the cell structures at about 5 and 1 8 mm from nozzle, the two sets follow a similar trend. 



Figure 9. Measurement of radial velocity, near the jet 
centerline slightly off axis, in a Mach 2 over-expanded jet. 
The red line plot shows an average of the CFD calculation 
at the measurement location. 


C. Axial velocity at two points near centerline 

Static measurements of velocity were performed using the type of interferogram shown in Fig. 4e. A region on 
the jet axis at 2.5 nozzle diameters upstream of the nozzle was selected for this experiment. The interferometer was 
set for a FSR of 5 GHz allowing a velocity measurement up to about 3500 m/sec without ambiguity. Four hundred 
interferograms containing information to compute three components of velocity were recorded at a laser pulse 
frequency of 20 Hz and analyzed as described earlier. The histogram of the measured axial velocity (V 3 ) at two 
locations, u 2 and u 3 situated 0.5 mm apart is presented in Fig. 10. 

The velocity fluctuations of about 35 m/sec are mainly due to the errors introduced by the on-line fitting 
method. Although the experimental noise of the facility was substantial, the simultaneous calculation of the 
reference reduced the errors introduced by vibrations detuning the interferometer to a minimal value. As before, the 


variation in the measured FSR (plotted around the 
full range of the instrument. Converted in velocity 
this uncertainty is about 25 m/sec. The 
uncertainty might possibly be caused by the errors 
in fitting the speckle-dominated reference. A 
potential enhancement would be to improve the 
reference light source and to use more elaborate 
fitting techniques of the Rayleigh-Brillouin 
spectra, which would take into account the 
interferometer transfer function and the intensity 
distribution of the combined input beam. Another 
way to improve the fitting technique precision is 
to increase the number of photon counts by 
increasing the laser energy or by increasing the 
photon efficiency of the camera. A more careful 
selection of the receiving optics, for an increased 
interferogram resolution and an increased number 
of photons per pixel received, would also help 
decrease these errors. 

The small difference in the measured velocity 
values at these two locations demonstrates the 
system’s capability for simultaneous 


mean as zero velocity) was found to be smaller than 1% of the 
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Figure 10. Instantaneous measurement of axial velocity 
at x = 25 mm on the axis of a Mach 2 over-expanded jet at 
P 0 = 5.5 atm and T 0 = 286.5 K. The two-point 
measurement of u 2 and u 3 is simultaneous. 
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measurement with high spatial resolution. In principle, the next outer fringe patterns can also be analyzed for 
velocity, but with reduced resolution. The average velocity at these two locations is about 5% higher than the value 
obtained by the CFD calculations, which is 422 m/sec. This discrepancy may be caused by uncertainty in the 
location of the measurement. 


IV. Conclusion 

A new technique using a single Fabry-Perot interferometer that enabled simultaneous, multiple-point 
measurements of up to three components of velocity was demonstrated. The velocity was determined by analyzing 
the interferogram containing spectral as well as spatial information. The mixing of the laser light and background 
scattered light with the Rayleigh scattered light in the same interferogram was used to reduce the velocity 
measurement uncertainty by simultaneous post processing. The adverse effects of the laser frequency drift, 
interferometer alignment errors, vibrations, and background light scattering were minimized by the use of this 
method. The laser energy at the measuring volume was increased about 4 times by pulse stretching, while keeping 
the energy per pulse under the limit of the gas breakdown. 

To approximate the Rayleigh-Brillouin spectrum, we used a combined function of one Gaussian and three 
Lorentzians in our fitting process to the measured data. For this system we showed that the fluctuations in the 
measured velocity (less than 25 m/sec) are mainly due to uncertainties in fitting the reference spectra. 

By determining slightly different velocity values at different locations in the flow, we demonstrated the system’s 
capability for simultaneous measurement with high spatial resolutions. 

With a few improvements, the technique is capable of obtaining accurate three-component velocity data at the 20 
Hz repetition rate required to improve CFD models. In combination with a dual-pump CARS technique, this method 
will provide single-pulse multi-property measurements for combustion research. 
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